Model for the hyperfine structure of electronically-excited KCs molecules 
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A model for determining the hyperfine structure of the excited electronic states of diatomic bialkali 
heteronuclear molecules is formulated from the atomic hyperfine interactions, and is applied to the 
case of bosonic 39 KCs and fermionic 4() KCs molecules. The hyperfine structure of the potential 
energy curves of the states correlated to the K(4s 2 Si/ 2 )+Cs( 6 p 2 A/ 2 , 3 / 2 ) dissociation limits is 
described in terms of different coupling schemes depending on the internuclear distance R. These 
results provide the first step in the calculation of the hyperfine structure of rovibrational levels of 
these excited molecular states in the perspective of the identification of efficient paths for creating 
ultracold ground-state KCs molecules. 

PACS numbers: 


I. INTRODUCTION 

The control of the evolution of atomic and molecu¬ 
lar systems at the single quantum level is an opportu¬ 
nity that is offered to fundamental physics by the amaz¬ 
ing developments in the research on ultracold quantum 
gases, i.e. where particles are moving with kinetic en¬ 
ergies E = ksT equivalent to temperatures T much 
smaller than 1 millikelvin. Among the most spectacular 
achievements are the observation of quantum degeneracy 
in ultracold atomic gases of bosons m and fermions 
[3]. Quantum degeneracy and even Bose-Einstein con¬ 
densation of weakly-bound ultracold molecules created 
by the association of pairs of atomic bosons or fermions 
has also been achieved EHTO], with the major drawback 
that dimer molecules are quite unstable against collisions 
as they contain a lot of internal energy, unless they are 
protected by the Pauli principle as in the case of fermion 
pairs. Major progress towards the formation of a quan¬ 
tum gases of ultracold molecules in their absolute ground 
state has been reported since 2008 Q3HH1], and the ob¬ 
servation of quantum degeneracy is now within reach. 
In most cases this research relies on two main steps: (i) 
the formation of weakly-bound molecules by association 
of a pair of ultracold atoms via a Feshbach resonance 
mm, and (ii) the transfer of the population from this 
so-called Feshbach state toward neighboring levels using 
avoided crossings m or radiofrequency transitions m , 
or toward the absolute ground level of the molecule us¬ 
ing a coherent optical process known as stimulated Ra¬ 
man adiabatic passage (STIRAP) [2TH25] . Alternative 
methods to obtain molecular degenerate gases are also 
under investigation starting from preexisting molecules, 
like evaporative cooling [23] , buffer gas cooling [25] 2B], 
Sisyphus cooling m or laser cooling |28H32] . 

All the above-mentioned methods require a precise 
control of the internal state of the molecule of interest. In 
particular the STIRAP technique has proven to be very 
versatile for controlling the internal state in all degrees 


of freedom (vibrational, rotational, and hyperfine) of ul¬ 
tracold diatomic molecules composed of identical atoms 
[TUI [T3] or of different alkali atoms [TT] [121 US]- This 
has allowed the creation of polar molecules, which are 
currently of great interest [33H36] . In brief, STIRAP is 
based on a so-called Lambda scheme of energy levels: the 
population is transferred from an initial level | i) toward 
a final level | g ) via two overlapping laser pulses involv¬ 
ing an intermediate level | e) that is never populated. 
For ultracold polar bialkali molecules, the scheme is im¬ 
plemented as follows: the | i ) level is a weakly-bound 
level of the molecular ground state manifold (usually re¬ 
ferred to as a Feshbach molecule) with mixed 3 S + and 
1 E + symmetries, the | g ) is the absolute ground level of 
the molecular ground state X 1 S + , and the | e) levels be¬ 
longs to an electronically excited state chosen such that 
it has noticeable dipole-allowed transition probabilities 
with both | i ) and | g). 

The knowledge of the tiniest properties of the quan¬ 
tum states involved in a STIRAP scheme - namely the 
hyperfine structure (hfs)- is essential to ensure the op¬ 
timal efficiency of the population transfer. The | i ) hfs 
is usually well-known for the alkali-atom pairs, at it is 
mostly determined by the hfs of the separated ground 
state 2 S' atoms involving only their electronic spin quan¬ 
tum number s = 1/2 and their nuclear spin quantum 
number i n . It has been largely used for the modeling of 
Fano-Feshbach resonances (FFR) in mixed alkali-metal 
atom pairs (i.e. LiNa [37], LiK [35], LiRb [33] @0], LiCs 
HUE2], NaK [33], NaRb [33], KRb [33], KCs [33], RbCs 
[47]). On the other hand very little is known on the hfs 
of ground-state heteronuclear alkali dimers in their low¬ 
est rovibrational levels (the | g ) state above), with a few 
remarkable exceptions [381153] . 

The hfs of the | e) levels is intrinsically more compli¬ 
cated as it involves one 2 P atom, i.e. one more non-zero 
angular momentum to be coupled. Its spectroscopic ob¬ 
servation could also be hindered by the natural width of 
the levels. Several studies have extended the long-range 
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analysis of the hfs of a pair of 2 S' ground-state alkali- 
metal atoms to the particular case where pure long-range 
potentials wells [SUES] are present in the electronic states 
dissociating to the lowest 2 S + 2 P asymptote of homonu- 
clear alkali dimers [B3H55] . Molecular spectroscopy has 
allowed the investigation of deeply-bound vibrational lev¬ 
els of excited states of NaRb [55 and NaK [S5HS5] . In 
the ultracold regime high-resolution spectroscopy has al¬ 
lowed observing and studying the hfs for low-lying vibra¬ 
tional levels of one of the first excited states of Rb 2 [63] 
and RbCs [61], which have been modeled with a simple 
effective Hamiltonian involving parameters fitted to the 
measurements. Last but not least, ab initio molecular 
hyperfine parameters have been computed for the lowest 
3 E + states of homonuclear alkali dimers [53 :. 

Despite all the achievements above, the complexity of 
the problem may have prevented up to now the com¬ 
plete description of the hfs of excited molecular states, 
as a prerequisite is the determination of complete po¬ 
tential energy curves (PECs) including hfs. In this 
paper, we extend the asymptotic model developed in 
Ref. [55] to all molecular states correlated to the two 
lowest 2 S + 2 P asymptotes of polar alkali-metal diatomic 
molecules. We illustrate it by deriving PECs includ¬ 
ing fine and hyperfine structure for all internuclear dis¬ 
tances R for the states of 39 KCs and 40 KCs correlated 
to K(4:S 2 Si/ 2 )+Cs(6p 2 Pj) (with j = 1/2, 3/2). We first 
recall the expression of the full Hamiltonian for the atom 
pair including spin-orbit and hyperfine couplings, and 
we detail the data used for the calculation for both the 
bosonic and the fermionic isotopologues. The variation 
of coupling regime along the PECs when R varies is dis¬ 
cussed in terms of ’’good” quantum numbers, shedding 
some light on the expected level structure depending on 
its binding energy. 

II. THE HAMILTONIAN FOR A PAIR OF K 
AND CS ATOMS 

Our goal is to compute Born-Oppenheimer (BO) PECs 
including fine and hyperfine structure. Thus the rotation 
of the molecule is not considered at that step. Within 
the standard BO approximation we write the electronic 
Hamiltonian of the system as 

H = H° + H f + H hf (1) 

where the electrostatic Hamiltonian H° containing the 
kinetic energy of the electrons and the sum of the poten¬ 
tial energies of charged particles is perturbed by the fine 
structure H f and the hyperfine structure H llf Hamilton 
operators. We assume the electrostatic problem solved 
so that the matrix of H° is diagonal and contains the 
P-dependent PECs in Hund’s case a (see Fig. [l] for the 
KCs curves correlated to K(4s 2 S)+Cs(6p 2 P)). 

For H f we only consider the molecular spin-orbit (SO) 
interaction, which is formally written as: 

H f = A(R) L.S (2) 


where L and S are the total electronic orbital momen¬ 
tum and spin. The P-varying coupling constant A(R) 
depends on the considered molecular states and is taken 
either from available spectroscopic data or from quantum 
chemistry calculations, as discussed below. 

Following Ref. [56l . we apply here the main approx¬ 
imation of our model: we express H llf as the sum 
of two atomic hyperfine Hamiltonian operators h hf (£:) 
(k = 1,2), which characterize the experimental hyperfine 
structure of the individual atoms. In other words, we 
assume that the related coupling constants do not vary 
with R , and we define them in order to reproduce the 
experimental hf splitting of the atoms. In the absence 
of published computations of the appropriate molecular 
constants -which are outside the scope of the present 
paper- this approximation is reasonable. For instance 
the variation of these coupling constants is expected to 
be limited to about ss15% in the case of RbSr ground 
state [66] ■ A similar trend is also found for the hfs of 
the 1 3 E+ and 1 3 E+ states of alkali-metal dimers .55] , 
Such an approximate Hamiltonian thus mimics the inter¬ 
action between the electrons and the nucleus within each 
atom usually expressed in terms of multipolar interac¬ 
tions, with odd magnetic and even electronic multipoles 
due to parity conservation inside the nucleus m-- 

h hf — h hf 4- L hf — L hf 4- L hf 4- L hf 4_ b hf f'Q'i 

11 — n MD T n EQ — n « + n si + n FC + n BQ W 

with the symbols £, s, and i referring to the atomic angu¬ 
lar momenta. The two main terms are: (i) the magnetic 
dipole interaction Hamiltonian which incorporates 

the interaction of the nuclear spin with the orbital an¬ 
gular momenta of the electrons h/f and with the spin of 
the electrons hk f , completed by the Fermi contact term 
hp f c induced by the electrons passing through the nuclei; 
(ii) the electric quadrupole term hpg, which describes 
the interaction between the electric quadrupole moment 
of the nucleus with the gradient of the electric field cre¬ 
ated by the electrons at the nucleus. In our model the 
interaction of the electronic and nuclear momenta of one 
atom with those of the other atom is neglected, so that 
we do not consider the corresponding multipolar terms. 

We apply the model to the first excited states of 
39 Ki 33 Cs and 40 K 133 Cs molecules, which tend to the 
K(4 2 5)-|-Cs(6 2 P) dissociation limit. Four electronic 
states are defined in the Hund’s case a notation: A 4 E + , 
c 3 E + , B 1 n and 6 3 n. Their potential energy curves 
(PECs) are presented on Fig. [I] They come either from 
the spectroscopic study of Refs. [55] [55], or from our own 
quantum chemical calculations based on effective core po¬ 
tentials (ECPs) and core polarization potentials (CPPs) 
implemented in a full configuration interaction approach 
(FCI) built on a large Gaussian basis set [7011751 . More 
details about the construction of the PECs can be found 
in m- Note that three out of four PECs dissociating 
to the next asymptote, K(4 2 P)+Cs(6 2 S) do not cross 
the PECs connected to the K(4 2 iS)+Cs( 6 2 P) lower limit. 
The (3) 4 E + and the P 1 n PECs cross each other but these 
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FIG. 1: (a) Hund’s case a potential energy curves of the 
excited electronic states correlating to the K(4 2 S)+Cs(6 2 P) 
(solid lines) and K(4 2 P)+Cs(6 2 S l ) (dashed lines) dissocia¬ 
tion limits. The B and c curves cross at R Bc = 7.65 a.u., 
the b and A curves at R bA = 9.55 a.u., and the b and c 
states at R bc = 12.84 a.u.; (b) the relevant P-dependent spin- 
orbit coupling terms associated to the curves dissociating to 
K(4 2 5 , )+Cs(6 2 P) (1 a.u. = a 0 = 0.0529177 mn). 


states have no direct spin-orbit interaction, and we will 
not include the former state in the present investigation. 

The chosen model for the interaction Hamiltonians 
above suggests to employ the well-known technique of 
linear combination of atomic orbitals to derive the ex¬ 
pressions of the coupling matrix elements of the spin- 
orbit and hyperfine interactions. We thus consider the 
basis \L, S, A, m$) where L and S are the total electronic 
angular momentum and spin quantum numbers, and A 
and ms the quantum numbers for their projection on the 
molecular axis: 

\L, S, A, ms) = * ( \hm 1 )^\l 2 m 2 )c s 

+ ( — 1 ) S l^l m l)K l^2W2)c s ) IS'ms) 1 ’ 2 . 

(4) 

The ket IS'ms) 1,2 is the coupled spin-function of the two 
valence electrons, and |/iTOi) and 1 / 27712 ) stand for the 
s and p atomic orbitals with l\ = 0, m\ = 0 and I 2 = 
1, m 2 = —1,0,1. For the present asymptotic limit we 
have L = 1 and A = m 2 , resulting in 12 basis functions 
\L,S,A,m s ). 


A. Spin-orbit coupling 

In the \L,S,A,ms) basis the H f operator couples 
states with the same A + ms = H value, where 
is the projection on the molecular axis of the total 
electronic angular momenta j e = L + S. As in the next 
step the hf interaction will couple most of the 12 states 
above, it is worthwhile to display the structure of H f 
operator in a compact and global way, as given in Eq. & 
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(5) 


For simplicity, this expression exhibits the location of the 
non-vanishing asymptotic (R -4 00 ) spin-orbit matrix el¬ 
ements, which have to be multiplied by the Cs(6p) spin- 
orbit splitting ( AE so /3 = 184.679 cm -1 ) [71] . The di- 
agonalization of this matrix yields the energies of atomic 
fine-structure levels. Each block is characterized by a 
value of Q (first line header in Eq. © , with the in¬ 
dex ± accounting for the reflection symmetry of = 0 
states) and involves Hund’s case a states (second line 
header), which are labeled as 2S+1 E™.£ 0 and 2S+1 n™ = f ±1 , 
with ms = 0, ±1. States with positive values of H are 
written in boldface. The degeneracy of the states with 
±fi 0 is not removed by the spin-orbit interaction. 

In order to obtain Hund’s case c PECs including spin- 
orbit, we have replaced the //-dependent spin-orbit ma¬ 
trix elements displayed in Fig. |ljb) to the atomic con¬ 
stant in Eq. ([ 5 ]). Note that all these functions indeed 
converge to AF so /3 = 184.679 cm” 1 at large distances. 
They are taken from various experimental and theoretical 
papers [BS], 75] [75], as explained in detail in Ref. [75]. The 
PECs resulting from the cliagonalization of the Hamilto¬ 
nian matrix including this molecular spin-orbit couplings 
are plotted in Fig. [2] To identify states with the same H 
value we label them with a short notation expressing their 
dissociation limit: 2(P 3 / 2 ), 1(Pl/ 2 ), 1(P 3 / 2 Q0> 1( P 3/2P), 
0+/”(Pi/ 2 ) and 0 +/ ”(P 3 / 2 ) notations. The two (P 3 / 2 ) 
states with Q = 1 are distinguished by the index a and 
(3 for the lowest and the upper one, respectively. 


B. Hyperfine interaction 

We build the hyperfine Hamiltonian upon the hyper¬ 
fine structure of the K and Cs atoms, which is recalled 
in Fig. [3] and in Table |IV| in the Appendix. The hyper¬ 
fine splitting of the 39,4 l *K(4 2 S 1 / 2 ) and the 133 Cs(6 2 Pi/ 2 ) 
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FIG. 2: Hund’s case c potential energy curves of ex¬ 
cited molecular states correlating to the dissociation limits 
K(4 2 5 , 1 / 2 )+Cs(6 2 Pi/2 ,3/2). The a and /3 indexes label the 
two fl = 1 states correlated to the upper dissociation limit. 
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state is determined only by the magnetic dipole term, the 
electric quadrupole and magnetic octupole terms being 
zero for j = 1/2. The latter terms give non-zero contri¬ 
bution only for the 133 Cs (6 2 P 3/ / 2 ) state. In the following 
we will omit the mass index of the Cs atom, for sim¬ 
plicity. The nuclear spins are 3/2, 4 and 7/2 for 39 K, 
40 K, and Cs, respectively. We set up the molecular basis 
by assembling the atomic states | ((sklk)jk,ik)fk,mf k ), 
where ik and fk denote the nuclear spin and the total 
angular momentum of the atom fc, respectively. We use 
the simplified notation 


|(/i,/2)/m/> 


\({{s\h)jih)h,({s 2 l2)j2i-i)f2)fm f ) 

{-l)~ h+h ~ mf ^/2f + 1 

m fl > m / 2 


X 

X 

X 


( fl h f \ 

V m /i m h ~ m f ) 

|((s2^2)j2,*2)/2,m/ 2 ), 


( 6 ) 


where (...) denotes 3j-symbols. The quantum numbers 
/ and m / are associated to the total angular momentum 
(without the molecular rotation) of the molecule f = fi + 
f *2 and its projection on the molecular axis. At infinite 
separation we obviously have 

1} /i) 

+ ^2 f (®2jj2, /2)>(7) 


E¥{ik,jk,fk) denoting the energy of the (■ h,jk,fk) hy- 
perfine level of the atom k. 

The matrices of H° and H f in this basis are ob¬ 
tained after a basis transformation between the asymp¬ 
totic molecular basis above and the \L,S,A,ms) ba¬ 
sis completed by the quantum numbers 1 1, mi) for the 
total nuclear spin of the molecule and its projection 


FIG. 3: Scheme of the asymptotic molecular hyperfine struc¬ 
ture of the bosonic * i9 KCs and fermionic 40 KCs molecules 
at the K(4 2 S' 1 / 2 )+Cs (6 2 Pi/ 2 j 3 / 2 ) dissociation limit. Energy 
spacings are not on scale, but the energy level order is mean¬ 
ingful. Exact energy splittings are reported in the Appendix. 


to/ = A + ms + mi on the molecular axis 


l(/i,/ 2 )/m/) = £ SLjijVjIfih 


j,I,rrij ,mi , 
S,L,ms ,rriL , 

,mi 2 


s^i3i rn j s-iL,m,L 

j,m,j,I,mi S,rris,L,mL 

ji h h 

h *2 /2 
J I f 

\L, S, A, m s ) |/to/), 



( 8 ) 


where X = \/2X + 1, the factors C are Clebsh-Gordan 
coefficients, and {...} are 9 j coefficients. We thus have 


1 = 2,...,5 for 39 KCs and I = 


f for 


40 


KCs. 


We have defined in Section IIA 12 | L, S, A, mg) states, 
while 32 (resp. 72) nuclear spin states | Imf) exist for 
39 KCs (resp. 40 KCs), resulting in a total of 384 (resp. 
864) |(/i 5 / 2 )/ to/ ) hyperfine states. As states with op¬ 
posite values of to/ are degenerate in energy in the ab¬ 
sence of external field, it is sufficient to count states 
with to/ > 0 (Tables [I] and [IT]). For instance, the num¬ 
ber of 39 KCs states with to/ = 0 is 48, so that there 
are (384-48)/2+48=216 different hyperfine states labeled 
with to/ > 0. As the hf interaction is implemented as 
a perturbation of Hund’s case c states, Tables [T] and [TT] 
also display the number of non-degenerate states for each 
fl > 0 symmetry. These numbers will be of great help 
to interpret the molecular PECs including hf interaction 
obtained in the next Section. 
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1 

13 

6 

2 

1 

1 

1 

0 

0 

0 

0 

5 

7 

1 

0 

0 

0 

0 

0 

0 

0 

1 

Total 

36 

36 

36 

36 

18 

18 

18 

18 

216 


TABLE I: The number of m/ > 0 hyperfine states in 39 KCs 
spread over the various |fi| symmetries. 
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Total 
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TABLE II: Same as Table Hj for 40 KCs. 


III. RESULTS 


We present our calculations for the hfs of the first ex¬ 
cited states of the bosonic 39 KCs and fermonic 40 KCs 
molecules, correlated to the K(4 2 S' 1 / 2 ,/i)+Cs(6 2 Pj 2 ,/ 2 ) 
dissociation limits displayed in Fig. [3} 

The long-range PECs resulting from the diagonaliza- 
tion of the total Hamiltonian (Eq. 0 ) in the molecu¬ 
lar basis of Eq. 0 are presented in Fig. [4] They are 
displayed down to R ~ 100 a.u. where the long-range 
van der Waals interaction is of the same magnitude than 
the hyperfine splittings. All PECs are attractive, as ex¬ 
pected from Fig. [5] Below 100 a.u. PECs correlated 
to different (/i,/ 2 ) asymptotes strongly recouple among 
each other, leaving / and to/ as the only remaining good 
quantum numbers. In this range the hfs is accurately 
calculated and could thus be of relevance for instance for 
photoassociation studies. Here we are mainly interested 
in the hyperfine structure of deeply-bound vibrational 
levels relevant for STIRAP, i.e. covering a range of lower 
distances. 

Below an internuclear distance of about 30 a.u., the 
Hund’s case c PECs are well separated in energy while 
their hyperfine structure is much smaller than their mu¬ 
tual spacing. Then it is convenient to plot the difference 
between each of these PECs with the related pure Hund’s 
case c ones. The global behavior of these difference PECs 
(or DPECs in the following) in the range on internuclear 
distances between 5 a.u. and 25 a.u. is displayed for both 
isotopologues in Figs. [7] and [8] (Q = 2,1 states) and in 
Figs. [5] and [8] (H = 0 states). For a given fl value, the 
DPECs exhibit similar amplitudes of the manifold, with 
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FIG. 4: Long-range potential energy curves (PECs) of 
KCs including hyperfine interaction, correlated to the 
K(4 2 5 1/2 ,/i)+Cs(6 2 P j2 ,/ 2 ). (a) and (c): the 216 PECs of 
39 KCs. (b) and (d): the 432 PECs of 40 KCs. Hyperfine dis¬ 
sociation limits are labeled with (/i,/ 2 ) 


zones of regularities and sharp transitions between them 
along the R interval, for both isotopologues. Most of 
such transitions result fr om the crossings between Hund’s 
case c PECs (see Table III) or Hund’s case a PECs (see 
Fig. 0 , which are then coupled by hf interaction so that 
the DPECs cannot be rigorously defined anymore. This 
pattern appears as more pronounced for the fl = 0 states, 
since the hyperfine splitting of these states are with 2 to 
3 orders of magnitude smaller than for D = 1,2 states. 
Note that the peak around 23 a.u. in the DPECs for 
0 + (P 3 / 2 ) and 0 _ (P 3/2 ) states, not reported in Table 


III 


is a consequence of the closeness (~ 0.068 cm -1 ) of these 
states at this location, so that their hyperfine structures 
overlap each other. 


In between these peaks, the DPECs may reflect the lo¬ 
cal Hund’s case a character of the hyperfine states, which 
is an important information for the future interpretation 
of the hyperfine structure of molecular levels. While our 
model does not explicitly treat the individual terms of 
the hfs Hamiltonian of Eqj3j one can discuss their influ- 
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1 (P3/2Q!) 

0 + (P3/2) 

0 (P3/2) 

2 ( P 3/2) 

12.82 

9.54 

12.81 

l(Pi/a) 


9.73 

12.66 


1(P3/2 0) 



12.80 

23.36 


TABLE III: Internuclear distances (in a.u.) where the KCs 
Hund’s case c PECs cross each other. Values in italics corre¬ 
sponds to the location of the sharp peaks in the corresponding 
difference PECs for (see text). 



FIG. 5: (Color online) Difference potential energy curves for 
= 0 states in 39 KCs, for every value of |m/| indicated by 
the color code. Each panel contains 18 curves (see Table [T|. 



m, 

1/2 
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5/2 

7/2 

9/2 

11/2 

13/2 

15/2 


FIG. 6: (Color online) Difference potential energy curves for 
fl = 0 states in 40 KCs, for every value of |m/| indicated by 
the color code. Each panel contains 36 curves (see Table [TT|). 


ence by formally rewriting in an effective manner for the 
molecule, so that the amplitude of the hfs manifolds are 
easily interpreted: 


H 


hf 

mol 


= H-lj + h 


hf 

SI 


frhf 
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rrhf 

n-EQ 


(9) 
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FIG. 7: (Color online) Difference potential energy curves for 
If = 2,1 states in 39 KCs, for every value of |m/| indicated by 
the color code. Each panel contains 36 curves (see Table [I]). 
When appropriate, the dominant Hund’s case a state is re¬ 
ported at the top of the panels for each interval between two 
crossings of Hund’s case a potential energy curves (indicated 
by solid vertical lines). The dashed lines at large R roughly 
locates the distance beyond which Hund’s case a states are 
strongly mixed by spin-orbit interaction. 



R(a.u.) 


FIG. 8: (Color online) Same as Fig. [ 7 ] for 40 KCs. Each panel 
contains 72 curves (see Table [TJ. 


where the symbols L , S and I now refer to the angular 
momenta of the molecule. 

The 2(-P 3 / 2 ) DPECs are solely built on the 6 3 II state, so 
that the resulting hfs is monotonous with R and is spread 
over about 0.8 GHz for 39 KCs (Fig. [tJi,) and 1.4 GHz for 
40 KCs (Fig. [ 8 ^,). The order of the [my curves is con¬ 
trasted between the two isotopologues: they are arranged 
in a somewhat disordered way in 40 KCs, while they ap¬ 
pears symmetrically around the my = 0 ones in 39 KCs, 
except for the largest |m/| values. This simple exam¬ 
ple illustrates the tight competition between the various 
terms of the hyperfine Hamiltonian, thus preventing from 
a uniform hfs pattern. In the 6 3 n state, the projection of 
the total electronic angular momentum as well as of the 
spin can not be zero. All of the three terms Hpj, H|j, 
and Hp f c of the magnetic dipole Hamiltonian do con- 
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tribute to the hfs, yielding a manifold with the largest 
amplitude of all the excited states. 

The variations of the 1 (P 7 ) DPECs display several 
regimes related to the local Hund’s case a character of 
the states. The |m/| curves are spread symmetrically 
around the to/ = 0 for the 1 (Pl/ 2 ) and l(P 3 / 2 ot) mani¬ 
folds in 39 KCs, and for the l(P %/2 ft) manifold in 40 KCs. 
In contrast, the DPECs are not structured according to 
|to/| values for the l(Pi/ 2 ) and l(p 3 / 2 <a) manifolds in 
40 KCs, and for the l(P 3 / 2 P) manifold in 39 KCs. The en¬ 
ergy amplitude spanned by the DPECs strongly depends 
on their local character, the largest one occurring when 
the states are mostly of c 3 E + type. The smallest hfs is 
found for the DPECs of P 4 n character for which the Hgj 
vanishes. A slightly larger hfs is visible for the DPECs of 
6 3 II character for which the projection of the electronic 
spin on the molecular axis vanishes in H|p The largest 
hfs is found for the DPECs of c 3 E + character as the three 
terms of the magnetic dipole interaction contribute to the 
interaction. 

The analysis of the structure of the DPECs for vari¬ 
ous ranges of distance reveals information on the main 
coupling regimes. In several cases, we observe that the 
DPECs are arranged in eight groups, corresponding to 
the possible values of the projection of the Cs nuclear 
spin nii Cs = —7/2,..., 7/2: the l(P 1 / 2 ) DPECs in the re¬ 
gion with 6 3 II character (Fig. 0 », Fig. [sjo) and the l(P 3 / 2 ) 
DPECs in the region with P 4 II character (Fig. [tJc, d 
and Fig. [8jc, d). Each of these eight groups exhibits a 
finer structure with DPECs arranged in either in four 
groups (Fig. [ 9 } for 39 KCs or nine groups (Fig. 101 for 
40 KCs molecule. This pattern correspond to the four 
values of rrii 3aK = —3/2,...,3/2 and to the nine values 
of mi iQK = —4,..., 4, respectively. This indicates that in 
the related R domains the nuclear spins of the two atoms 
do not couple each other. In 39 KCs, the l(Pi/ 2 ) DPECs 
between 13.8 a.u. and 16 a.u. and the l(P 3/2 ct) DPECS 
at 14.5 a.u., mostly of c 3 E + character are grouped in 11 
subsets corresponding to the projections of the total nu¬ 
clear spin To/ = —5,..., 5. These subsets are equidistantly 
located around the to/ = 0 one, revealing the dominant 
first order magnetic dipole interactions yielding energies 
proportional to to/ m- 



FIG. 9: Zooms on the it-intervals for which the 39 KCs 1(P/) 
DPECs of Fig. [7] are grouped according to the projection of 
the nuclear spin of the cesium and potassium atoms. The 
figure displays the four curves rm 39K = —3/2,..., 3/2 associ¬ 
ated to the 4th pack (counting from the top) in case of the 
39 K 133 Cs molecule. 

The coupling of the nuclear spins appears in other 
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FIG. 10: Same as Fig.^for the 40 KCs 1(P 7 ) DPECs of Fig. |] 
with nine curves rm i0K = —4, ...,4. 


cases. In 39 KCs, the 1 T \/ 2 DPECs between 13.8 a.u. 
and 16 a.u. and the l(P 3 / 2 o:) DPECS at 14.5 a.u., 
mostly of c 3 E + character are grouped in 11 subsets cor¬ 
responding to the projections of the total nuclear spin 
to/ = —5,..., 5. These subsets are equidistantly located 
around the to/ = 0 one, revealing the dominant first 
order magnetic dipole interactions yielding energies pro¬ 
portional to TO/ m- Furthermore, Fig. [ 8 Jd illustrates the 
coupling of nuclear spins for R > 16 a.u. in the l(Pi/ 2 ) 
DPECS of 40 KCs with sixteen groups corresponding to 
the possible projections of the total nuclear spin 15/2. 
This is easily understood as unlike 39 KCs, the hfs of the 
39 K(4s) is larger than the Cs( 6 pi/ 2 j 3 / 2 ) ones, so that the 
grouping of DPECs is different. 

In contrast with the cases above, the hyperfine energy 
splittings for the ST = 0 DPECs have a magnitude in the 
MHz range for both isotopologues, which can be inter¬ 
preted with Eq. [9] The first-order magnetic dipole inter¬ 
action vanishes, but not the second-order magnetic dipole 
and higher-order multipole interactions, which give small 
corrections to the energy EZj. 



R(au) 


FIG. 11: Zooms of the 39 KCs 0 + (P 3/2 ) DPECs of Fig. [HJr for 
ranges of R in between the peaks marking the transition from 
one Hund’s case a character to the other. 
































FIG. 12: Zooms of the 39 KCs CT(P 3/2 ) DPECs of Fig. [Hj:l for 
ranges of R in between the peaks marking the transition from 
one Hund’s case a character to the other. 
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FIG. 13: Zooms of the 40 ICCs 0 + (P 3/2 ) DPECs of Fig. [e}> for 
ranges of R in between the peaks marking the transition from 
one Hund’s case a character to the other. 



FIG. 14: Zooms of the 40 KCs CT(P 3/2 ) DPECs of Fig. [eji for 
ranges of R in between the peaks marking the transitionfrom 
one Hund’s case a character to the other. 


IV. DISCUSSION: TOWARDS A FULL MODEL 
OF THE FORMATION OF ULTRACOLD KCS 
MOLECULES 


We have extended the asymptotic model of Ref. (HU 
down to small internuclear distances in order to deter¬ 
mine the potential energy curves of electronically excited 
molecular states of 39 KCs and 40 KCs including hyperfine 
structure. The hyperfine splittings of the 0 = 1,2 molec¬ 
ular states have an amplitude in the GHz range, and 
typically three orders of magnitude less for the 0 = 0 
states, as expected from Ref. m- These calculations 
must be complemented by introducing molecular rotation 
and thus evaluate the hyperfine structure of the rovibra- 


tional levels of the molecules, which will be the purpose of 
a forthcoming paper. Nevertheless, the present calcula¬ 
tions are already relevant for various aspects of ultracold 
molecule research. 

First, the implementation of optical shielding in ul¬ 
tracold gases, i.e. the suppression of inelastic collisions 
between ultracold atoms, relies on the modification of 
long-range interactions between ground-state and excited 
atoms by laser light m- The experimental proof-of- 
principle was first achieved in trap-loss experiments in¬ 
duced by inelastic collisions of cold alkali atoms [75U5T| . 
The suppression of ionizing (Penning) collisions was ob¬ 
served with Xe [52] and Kr atoms [521, as well as the sup¬ 
pression of reactive collisions like photoassociative colli¬ 
sions of sodium atoms |84j . These experiments have been 
achieved for gases with temperature in the 10-100 pK 
range, too large to allow for selecting a single repulsive 
channel so that the suppression effect is partial. A full 
description of this process at even lower temperatures 
must include the hyperfine structure of the related elec¬ 
tronic states as determined in the present work. Here all 
the electronic states are attractive at large distances, but 
a few electronic states correlating to the K(4p)+Cs(6s) 
dissociation limit are indeed repulsive (see Fig. |T|) and 
thus suitable for optical shielding. 

Second, the most advanced achievements toward the 
creation of dense gases of cold molecules in their ground 
state rely on two main steps First 

a pair of ultracold atoms is associated into a loosely- 
bound vibrational level of the ground molecular state 
manifold using magnetic or laser field. Then a con¬ 
trolled population transfer from this well-defined level 
to the absolute ground state level can be coherently per¬ 
formed by the method of stimulated rapid adiabatic pas¬ 
sage (STIRAP), involving an intermediate excited molec¬ 
ular state. The theoretical modeling of such a multistep 
process without taking in account the hyperfine struc¬ 
ture of the molecular levels yields an efficient guiding 
procedure for experimental work [731 EH] However, the 
optimal efficiency of the transfer imposes the control of 
the population at the single quantum state level, namely 
including the hyperfine structure. The present results de¬ 
liver some hints for choosing a convenient intermediate 
state for STIRAP. This is exemplified by Fig. [15] show¬ 
ing the hyperfine structure of all the 39 KCs and 40 KCs 
PECs dissociating into K(4s 2 5' 1 / 2 )+Cs(6p 2 P 1 / 2 ,3/2) at 
Ri = 10.5 a.u.. This particular distance is chosen after 
the theoretical study of Ref. m, which proposes various 
options for achieving STIRAP in KCs involving interme¬ 
diate levels with radial wavefunctions located around Ri. 
The corresponding rotational energy estimated from the 
rotational constant B = h 2 /(2pRf) amounts for about 
500 MHz. The H = 0 vibrational levels will exhibit a 
well-defined rotational structure with each sublevel com¬ 
posed of a tiny hyperfine manifold, which may be te¬ 
dious to resolve experimentally. On the other hand, the 

= 1 vibrational levels will be characterized by a com¬ 
plex structure resulting from the competition of the hy- 
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perfine interaction and the molecular rotation of compa¬ 
rable magnitude. As already demonstrated experimen¬ 
tally, such levels can be addressed individually, but their 
characterization in terms of quantum numbers is more 
involved. For an individual level, it is possible to set up 
an effective Hamiltonian whose parameters are fitted on 
the recorded spectrum m In our forthcoming work, we 
will model the structure of all levels starting from the 
atomic hyperfine parameters, and we will adjust them in 
order to match the observations. 
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Appendix 

We present in Table |IV| the values of the hyperfine 
energy splittings of the 3y K and 40 K 4s level, and of the 
Cs 6 p 1 2 3 Pi/ 2 , 3/2 levels in terms of the related molecular 
dissociation limits. 


K(/s)+Cs(6p) 

(IaAcs) 

K (4sj K )+Cs(6p jc J 

(j A, j Cs ) 

E/(cm 4 5 6 7 8 9 ) 

39 K (/ Sjk )+ 133 
Eh/(cm 4 ) 

Cs (6Pics) 

(fcsjfx) 

40 K(/^) + 

Eh/ (cm 4 ) 

133 Cs(f>p, c J 

(fcs,fKr) 




184.7498 

(5,2) 

184.712 

(5,7/2) 




184.7415 

(4,2) 

184.704 

(4,7/2) 




184.7345 

(3,2) 

184.697 

(3,7/2) 


(1/2, 3/2) 

184.7351 

184.7341 

(5,1) 

184.692 

(2,7/2) 




184.7294 

( 2 , 2 ) 

184.67 

(5,9/2) 




184.7258 

(4,1) 

184.661 

(4,9/2) 




184.7191 

(3,1) 

184.654 

(3,9/2) 

( 0 , 1 ) 



184.7141 

( 2 , 1 ) 

184.649 

(2,9/2) 




-369.4463 

(4/2) 

-369.319 

(4,7/2) 


( 1 / 2 , 1 / 2 ) 

-369.4696 

-369.4629 

(4,1) 

-369.358 

(3,7/2) 




-369.4863 

(3,2) 

-369.362 

(4,9/2) 




-369.5030 

(3/T) 

-369.401 

(3,9/2) 


TABLE IV: Asymptotic molecular hyperfine ener¬ 
gies of the 39 K 133 Cs and 40 K 133 Cs molecules for the 
K(4s 2 Si/ 2 )+Cs( 6 p 2 Pi/ 2 , 3 / 2 ) dissociation limits, see fig. [ 3 ] 
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FIG. 15: Difference potential energies of the Hund’s case c 
states at R = 10.5 a.u. of the bosonic 39 KCs (left panels) and 
fermionic 40 KCs (right panels) molecules. The lower panels 
display zooms for the Q = 0 states. Both the relevant Hund’s 
case c states and the dominant Hund’s case a character are in¬ 
dicated. For Hund’s case c notation in the upper line we have 
used a simplified one giving only the value of D. From column 
1 to column 8 the order of the states is l(-P 3/2 ck), 0 “(P 3 / 2 ), 
2 (f3/2) 5 1(T’i/ 2 ), 0 (P 1 / 2 ), 0+(P 1/2 ), l(P 3 /2 P), and 0 + (P 3 / 2 ). 





























